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GENERATING  SHIP  MOTION  TIME  HISTORIES 
By 

P  Crossland 


Summary 

PAT-86  is  a  suite  of  seakeeping  computer  programs  which  may  be  used  for 
predicting  motion  spectra  of  ships  in  waves.  This  Technical  Memorandum 
describes  a  method  of  generating  time  histories  from  these  motion  spectra. 
A  computer  program  which  generates  time  histories  has  been  written  and  is 
listed  in  the  Appendix. 

Comparisons  of  the  motion  spectra  from  PAT-86  are  made  with  the  spectra  of 
the  generated  time  histories.  These  comparisons  are  shown  to  be 
satisfactory. 
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Notation 


6w  Frequency  interval 

Sio„  Encounter  frequency  interval 

£ 3  Total  motion  phase 

£ r  Wave  phase 

£ z  Heave  phase 

g  Acceleration  due  to  gravity 

fi  Heading 

N  Number  of  frequency  intervals 

sno  Total  motion  amplitude 

S,  Pitch  spectrum 

Sf  Wave  spectrum 

S2  Total  motion  spectrum 

s(t)  Total  motion  time  history 

Wave  amplitude 
5"(t)  Wave  time  history 

t  time 

Th  Signal  repeat  period 

9{ t)  Pitch  time  history 

U  Speed 

w  Wave  frequency 

w.  Encounter  wave  frequency 

xB1  Distance  from  LCG  to  point 

of  total  motion  calculation 

zno  Heave  amplitude 

z(t)  Heave  time  history 


[ Rad/s ] 
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[Rad] 
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GENERATING  SHIP  MOTION  TIME  HISTORIES 


by  P  Crossland 

1.  OBJECTIVES 

There  has  been  a  recent  demand  for  time  histories  of  ship  motions  to  be 
calculated  for  a  given  ship  in  a  given  sea  condition.  The  PAT- 86  suite  of 
ship  motion  computer  programs  calculates  motion  spectra  rather  than  the 
instantaneous  values  or  time  history. 

This  technical  memorandum  outlines  a  method  of  generating  time  histories 
from  the  five  principle  ship  motion  spectra  and  from  the  total  motion 
spectra.  The  problems  associated  with  generating  multiple  time  histories 
are  related  to  the  choice  of  phases .  Random  phases  must  be  chosen  in  such 
a  way  that  the  correct  correlation  between  the  motions  is  maintained. 

Examples  of  time  histories  are  shown  and,  as  one  check,  the  motion  time 
history  at  a  point  is  calculated  from  the  heave  and  pitch  motion  and 
compared  with  that  generated  from  a  total  motion  spectrum. 

2.  INTRODUCTION 

As  the  ship  passes  through  the  waves  it  responds  in  six  rigid  body  modes, 
surge,  sway,  heave,  roll,  pitch  and  yaw.  The  sign  convention  and  units  are 
given  in  Reference  1.  These  rigid  body  motions  are  superimposed  on  one 
another  leading  to  a  compound  motion  at  any  point  on  the  ship. 

These  rigid  body  motions  and  compound  motions  are  calculated  by  the  PAT-86 
suite  of  computer  programs  at  ARE  (Haslar) .  The  suite,  first  developed 
around  1980,  was  based  on  state-of-the-art  theory.  The  prediction  of  ship 
motion  is  calculated  from  frequency  domain  dynamic  analysis  of  the  ship  in 
waves  and  the  process  is  divided  into  three  stages . 

a.  Calculation  of  motion  in  regular  (sinusoidal)  waves  over  a  range 
of  frequencies  of  interest.  This  calculation  uses  strip  theory. 

b.  Combination  of  these  regular  wave  responses,  or  transfer 
functions,  with  a  wave  spectrum  to  give  a  response  spectrum. 

c.  Integration  of  the  response  spectrum  to  give  variance  of  motion 
displacement.  Calculation  of  the  second  and  fourth  moments  will  give 
variances  of  velocity  and  acceleration  respectively. 

A  full  description  of  the  theory  is  given  in  Reference  2.  This  TM  outlines 
a  method  of  converting  the  response  spectrum  into  time  histories  of  that 
response.  The  method  used  is  summation  of  sine  waves  with  randomly 
generated  phases. 

3.  PAT-86  OUTPUT 

The  output  from  the  PAT-86  suite  contains  the  regular  wave  calculation 
(transfer  function  and  phase)  and  the  irregular  wave  calculation  of  the  5 
rigid  body  motions  and  total  motions  at  specified  points.  (Surge  is 
neglected) .  All  this  information  is  required  if  time  histories  are  to  be 
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generated.  So  it  must  be  extracted  from  the  main  PAT-86  output  file  by 
some  suitable  means . 

4.  WAVES 

A  wave  spectrum  is  input  into  PAT-86  and  the  resulting  motion  spectra  are 
calculated.  The  wave  spectrum  is  in  the  wave  frequency  domain  and  the 
motion  spectra  are  in  the  encounter  frequency  domain. 


The  first  step  in  generating  time  histories  is  to  change  the  wave 
spectrum  to  a  spectrum  in  the  wave  encounter  frequency  domain.  This 
process  is  known  as  warping. 

The  wave  encounter  frequency  becomes : 

„  -  „  _  u2uccs( u)  ... 


The  spectral  ordinate  becomes: 


S{(coe) 


Sc(<a) 

2  a)  Ucos  ( |i ) 


The  Power  Spectra  output  from  PAT-86  need  to  be  converted  to  amplitude 
spectra  using  the  following  formula  for  the  wave  amplitude  associated 
with  the  n0*  frequency: 


(<*),) 


where 


c .  Generation 

The  time  history  is  generated  by  adding  together  a  large  number  of 
component  sine  waves  whose  amplitude  is  given  by  the  amplitude 
spectrum.  Thus,  the  time  history  is  given  by: 

C  ( t)  =  E  cos  (uet  +  ec)  ( 


for  0  <  t  <  T„ 


2 


(6) 


where  TH  is  the  signal  repeat  period  given  by 

T  _  2n 
1h  " 


and  is  a  randomly  chosen  phase  lying  between  -n  and  n  (different  for 
each  frequency) . 

5.  RIGID  BODY  MOTIONS 

PAT-86  outputs  the  rigid  body  motions  with  phase  relative  to  a  wave  peak  at 
the  CG  of  the  ship.  This  calculation  is  carried  out  in  regular  waves.  The 
method  of  synthesising  time  histories  of  the  rigid  body  motions  is  similar 
to  that  described  in  Section  4.  However  the  regular  wave  phases  are  added 
to  the  random  phases  generated  in  4. 

Thus,  taking  heave  as  an  example. 

The  heave  amplitude  associated  with  the  nth  frequency  is  given  by 

Zno  =  v/2S2(ue)  6o>„ 


and  the  time  history  is  given  by 

w 

z(t)  =  2  z„0cos(wet  *  (e{  +  £j)) 

n-1 


for  0  <  t  <  Th 

c z  is  the  heave  phase  taken  from  the  regular  wave  response  (transfer 
function).  The  other  4  motions  (sway,  roll,  pitch  and  yaw)  are  dealt  with 
in  a  similar  way. 

6.  TOTAL  MOTIONS 

6.1.  Generation  using  PAT-86  Total  Motion  Spectra 

PAT-86  outputs  the  total  motions  (absolute  vertical , lateral  and  relative 
vertical)  at  specified  points.  Again  the  phases  cz  are  found  in  regular 
waves  and  these  must  be  combined  with  the  random  phases  c f  found  in  4,  to 
generate  a  time  history  from  the  calculated  motion  spectrum.  The  motion 
amplitude  spectrum  is  found  from 

s„0  =  6<o 7  (9) 


and  the  resulting  time  history  is  given  by 

N 

s{c)  =  2  sn0cos(ti> et  +  (e{  +  «,)) 

n-1 


(10) 
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6.2.  Generation  using  Rigid  Body  Motion  Time  Histories 


An  alternative  method  of  calculating  total  motion  time  histories  is  now 
described.  Generate  the  time  histories  of  the  five  rigid  body  motions 
using  the  method  outlined  above  in  Section  5.  Then,  combine  the  time 
histories  to  obtain  the  total  motions  at  a  point. 

Take  the  simple  case  of  vertical  motion  at  a  point  xB1  metres  forward  of 
the  LCG  and  on  the  ship  centreline.  The  time  history  of  vertical  motion  is 
given  by 

Sit)  =  zit)  -  XB1d  it)  (11) 

where  the  pitch  angle,  9,  is  in  radians. 

7.  EXAMPLES 

A  ship  motion  calculation  was  carried  out  using  the  PAT-86  suite. 

The  calculation  considered  only  the  vertical  plane  motions  of  a  ship 
travelling  in  head  seas. 

The  program  calculated  the  absolute  vertical  motion  at  the  FP  on  the 
centreline  of  the  ship.  Time  histories  of  heave  and  pitch  were  calculated 
and  the  resulting  total  motion  time  history  was  determined  using 
equation  (11).  As  a  check,  a  total  motion  time  history  was  generated  from 
the  motion  spectrum  using  equation  (10).  Comparisons  were  made  of  both  the 
total  motion  time  histories  and  their  spectra. 

A  program  was  written  to  generate  the  time  histories  and  carry  out  the 
comparisons.  The  program  is  shown  in  Appendix  A. 

As  a  further  study,  correlation  coefficients  of  the  heave  and  pitch  time 
histories  were  calculated  in  accordance  with  the  method  described  in 
Reference  3. 

8.  RESULTS 

Figure  1  shows  the  wave, heave  and  pitch  spectra,  found  in  PAT-86,  from 
which  the  time  histories  were  generated.  Also  shown  are  the  spectra 
calculated  from  the  generated  time  histories  to  confirm  that  they  are 
correct.  The  small  differences  between  the  two  spectra  are  probably  due  to 
the  short  signal  repeat  period  used.  Time  histories  from  the  spectra  are 
shown  in  Figure  2.  A  50  second  section  of  each  time  history  is  shown,  but 
about  1200  seconds  was  generated  in  all.  Figure  3  shows  the  variation  of 
Heave-Pitch  correlation  coefficient  with  heading.  The  curve  calculated 
from  the  generated  time  histories  is  compared  with  the  correlation  curve 
found  from  the  seakeeping  trial  (Reference  3).  The  trends  shown  in  both 
curves  are  similar,  however  the  actual  values  are  quiue  different.  This  is 
not  particular  surprising  since  two  different  ships  are  being  compared,  but 
it  illustrates  that  the  time  histories  are  correctly  correlated.  This 
would  suggest  that  the  phase  information  has  been  chosen  correctly. 

Figure  4  shows  total  motion  spectra  at  the  FP.  The  total  motion  spectrum 
as  found  in  PAT-86  is  compared  with  the  total  motion  spectrum  calculated 
from  the  heave  and  pitch  time  histories.  The  small  differences  in  the  two 
spectra  may  again  be  attributed  to  the  complete  time  history  being  two 
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shore.  Finally,  the  total  motion  time  history  generated  from  the  total 
motion  spectrum  is  shown  in  Figure  5.  The  total  motion  time  history  found 
from  heave  and  pitch  is  also  shown  but  it  is  actually  coincident  with  the 
first  time  history. 

9.  CONCLUSION 

This  Technical  Memorandum  describes  a  method  of  generating  time  histories. 
The  comparisons  shown  indicate  that  the  choice  of  phasing  the  random  signal 
is  correct. 
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Appendix  A 


FORTRAN  PROGR^ 


PROGRAM  HISTORY 

* 

*  Program  to  generate  time  history  given  a  spectrum. 

*  The  phase  information  is  taken  from  a  uniform  random 

*  variable  between  -PI  AND  PI 


*  Written  by 

*  P  Crossland  04-0ct-1990 

* 

* 


REALMS  X.G05CAF 

REAL  XG( 3200) , A( 50) , B( 50) , SEC( 50) . W( 50) , S( 50) ,ZETA(50) 

*  , E( 50) ,HEAD(2> ,WE(50) ,SE(50) ,XGS(3200) ,XGC(3200) 

*  ,  TME( 3200 ) , DWE( 50 )  SC ( 50) , EZ( 50) , ETHETA( 50 ) 

*  ,  Sz.(  50 )  ,  STHETA(  50)  ,  ES(  50 ) 

*  ,SS(50),C(3200), TOT ( 50 ) . HEAV( 50 ) , PITC( 50 ) , Z( 3200 ) 

*  , THETA( 3200) ,FP( 3200) 

COMMON/DATA/NTIME , WE , TIME . DWE , ANS 

* 

*  "''ta  statement  containing  necessary  variables  which  must  be 

*  changed  according  to  the  PAT-86  output. 

* 


DATA  OMEGA I , OMEGAF , DW , U  GRAV  /0 . 05 , 2 . 50 , 0 . 05 , 10 . 296 , 9 . 81/ 
DATA  SUM4/5*0 . 0/ 

DATA  HEAD  /45.0, 180.0/ 

* 

*  The  information  contained  in  the  above  data  statements  can  be 

*  obtained  from  the  output  file  from  PAT- 86 

* 

* 

*  Open  input  file  WAVE.O  containing  a  maximum  of  50  spectral 

*  ordinates  with  frequency  as  the  independent  variable. 

* 


OPEN ( UNIT=C , FILE- ' WAVE . 0 ' ) 


* 

*  Open 

rest  of  the  f: 

* 

REG .  0 

* 

* 

IRREG. 0 

* 

TREG . 0 

★ 

TIRREG. 0 

contains  the  phase 
pitch  and  heave 
contains  the  heave 
contains  the  phase 
contains  the  total 


information  for 

and  pitch  spectra 
information  for  motions 
motion  spectrum 


calculated  spectra 


OPEN(  (JNIT-7  ,  FILE” '  REG  .O') 

OPEN( UNIT-9,  FILE- 'IRREG. O' ) 

0PEN( UNIT-10,  FILE- ' TREG.O' ) 

OPEN( UNIT-11,  FILE- 'TIRREG. O' ) 

* 

*  Open  outp"t  files  SPECTRUM  and  TIME,  the 

*  and  the  calculated  time  histories 
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* 

0PEN(UNIT=8, FILE-' TIME ' ) 

OPEN ( UNIT-6 , FILE— ' SPECTRUM' ) 

* 

*  Initialise  variables  and  declare  constants. 

* 

PRINT*, 'What  is  the  heading  in  degrees' 

READ(* , * )  HEAD( 2 ) 

PI  =  4 . 0*ATAN( 1.0) 

NPOINTS  =  0 
IWT  -  0 
IFAIL  -  0 

* 

* 

*  Read  in  the  frequencies  and  spectral  ordinates  and  phase 

*  information 

*  D1 . D2 . DJ . D4 . ID  are  insignificant  information. 

* 

DO  10  1=1,50 

READ( 5,*, END-11)  W(I),S(I) 

READ( 7  .  * , END-11 )  D1 , D2 . D3 , EZ( I ) . D4 , ETHETA( I ) 

READ( 9,*, END-11)  ID , D2 , SZ( I ) , STHETA( I ) 

READ( 10.*, END-11)  DL . D2 , D3 . ES( I ) 

READ( 11,*, END-11)  ID , D2 , SS( I ) , D3 , D4 

* 

*  Change  phases  into  radians  and  count  the  number  of  frequency 

*  steps. 

* 

EZ( I )  -  EZ( I )*PI/180 . 0 
ES ( I )  -  ES( I )*PI/180 . 0 
ETHETA(I)  =  ETHETA( I )*PI/180 . 0 
NPOINTS  -  NPOINTS  +  1 

10  CONTINUE 

11  CONTINUE 

* 

*  Warp  the  wave  spectrum  to  encounter  frequencies 

*  Change  Heading  input  to  radians 

* 

HEAD(l)  -  HEAD(1)/180.0*PI 
HEAD( 2 )  -  HEAD ( 2 )/180 . 0*PI 
DO  100  I  -  1, NPOINTS 

WE( I )  -  W(I)  -  W(I)*W(I)*U*C0S(HEAD(2))/GRAV 
SE( I )  -  S ( I ) *GRAV/ ( GRAV - 2 . 0*W( I )*U*C0S(HEAD( 2 ) ) ) 

100  CONTINUE 

Cl  -  2 . 0*U*C03 ( HEAD( 2 ) )/GRAV 
C2  -  -  2 . 0*U/GRAV 

* 

*  Calculate  the  amplitude  spectrum  for  the  required  range 

*  of  frequencies. 

* 

DO  12  I  -  1, NPOINTS 
DWE(I)  -  a.0-Cl*W(I))*DW 
TOT(I)  -  SQRT( 2 . 0*SS( I )*DWE( I ) ) 

HEAV(I)  -  SQRT( 2 . 0*SZ( I )*DWE( I ) ) 
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PITC(I)  =  SQRT(2 . 0*STHETA( I )*DWE( I ) ) 

ZETA(I)  =  SQRT  < 2 . 0*SE( I )*DWE( I ) ) 

12  CONTINUE 

* 

*  Calculate  the  repeat  period  of  the  signal  to  be  generated 

* 

TH  =  2 . 0*PI/DW 


* 

*  Generate  the  random  phases  using  NAG  libraries  and  infer  the 

*  phases  for  heave .pitch  and  total  motion 

* 

CALL  G05CBF(0) 

DO  13  I  =  1 , NPOINTS 
X  =  G05CAF(X) 

E( I )  -  (FLOAT(X) -0 . 5)*PI*2 . 0 
EZ( I )  =  E ( I )  +  EZ( I ) 

ETHETA(I)  -  E ( I )  +  ETHETA(I) 

ES ( I )  =  E( I )  +  ES ( I ) 

13  CONTINUE 

* 

*  Generate  the  time  history  for  TH  seconds  with  an 

*  time  interval  of  0.25  seconds 

* 

NT I ME  -  4 . 0*NINT(TH) 

DO  14  I  -  1,  NTIME+1 
T  -  FLOAT( (I-l)/4.0) 

SUM  -  0.0 

SUM1  =0.0 

SUM2  =0.0 

SUM3  =0.0 

DO  15  J  =  1, NPOINTS 

SUM  -  ZETA ( J ) *C0S ( WE ( J ) *T  +  E(J))  +  SUM 
SUM1  =  HEAV( J )*C0S(WE( J)*T  +  EZ(J))  +  SUM1 
SUM2  -  PITC( J)*C0S(WE( J)*T  +  ETHETA(J))  +  SUM2 
SUM3  -  TOT ( J ) *C0S ( WE ( J ) *T  +  ES(J))  +  SUM3 
15  CONTINUE 

XG( I )  -  SUM 
Z(I)  =  SUM1 
THETA ( I )  =  SUM2 
FP( I )  =  SUM3 

* 

*  Calculate  the  total  motion  at  the  FP  from  the  assumed  time 

*  histories 

C(I)  -  Z(I)  -  64 . 20*THETA( I )*PI/180 . 0 

* 

*  Ouput  the  time  histories 

* 

WRITE! 8. 1001)  T,C!I),FP!1).Z!I),THETA(I),XG(I) 

14  CONTINUE 

1001  FORMAT! F8 . 4 , 5F12 . 6 ) 
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*  Calculate  the  spectrum  of  the  generated  time  history  to  check 

* 

*  First  calculate  the  spectral  ordinates 

* 

DO  200  I  =  1 , NPOINTS 
CALL  F0URIER( Z , I ) 

HEAV(I)  =  ANS 

CALL  FOURIER (THETA, I ) 

PITC(I)  =  ANS 
CALL  FOURIER(FP, I) 

TOT ( I )  =  ANS 
CALL  FOURIER ( C , I ) 

SC(I)  -  ANS 
CALL  FOURIER (XG , I ) 

S(I)  =  ANS 

* 

*  Ouput  the  spectra 

* 

WRITE( 6 , 1000)  WE( I ) , HEAV( I ) , PITC( I ) , TOT( I ) , S( I ) , SE( I ) 

200  CONTINUE 
1000  FORMAT ( 6F12 . 5 ) 

STOP 

END 

SUBROUTINE  FOURIER( ZZ . K) 

REAL  ZZ( 3200) , WE( 50) , TIME( 3200 ) ,XGC(3200) ,XGS(3200) ,DWE(50) 
COMMON/DATA/NTIME , WE , TIME , DWE , ANS 
SUM1  -  0.0 
SUM2  -  0.0 
DO  210  J  =1 , NTIME 
TIME(J)  =  FLOAT(J-1)/4.0 
XGS(J)  -  ZZ( J)*SIN(WE(K)*TIME( J) ) 

SUM1  -  SUM1  +  XGS(J) 

XGC(J)  -  ZZ(J)*COS(WE(K)*TIME(J)) 

SUM2  -  SUM2  +  XGC(J) 

210  CONTINUE 

A  =  SUM2*2.0/FLOAT( NTIME) 

B  -  SUM1*2 .0/FL0AT( NTIME) 

* 

*  Calculate  the  spectral  ordinates 

* 

ANS  -  (A*A+B*B)/2.0/DWE(K) 

RETURN 

END 
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